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> 1 INTRODUCTION 



ABSTRACT 

We report the first IR proper motion measurements of the Herbig-Haro objects in 
the Orion Molecular Cloud-One using a four-year time baseline. The [Fell] emitting 
bullets are moving of order 0.08 arcsec per year, or at about 170 knis^^. The direction 
of motion is similar to that inferred from their morphology. The proper motions of 
H2 emitting wakes behind the [Fell] bullets, and of newly found H2 bullets, are also 
measured. H2 bullets have smaller proper motion than [Fell] bullets, while H2 wakes 
with leading [Fell] bullets appear to move at similar speeds to their associated bullets. 
A few instances of variability in the emission can be attributed to dense, stationary 
clumps in the ambient cloud being overrun, setting up a reverse-oriented bullet. Dif- 
ferential motion between [Fell] bullets and their trailing H2 wakes is not observed, 
suggesting that these are not separating, and also that they have reached a steady- 
state configuration over at least 100 years. The most distant bullets have, on average, 
larger proper motions, but are not consistent with free expansion. Nevertheless an 
impulsive, or short-lived 1,000 years) duration for their origin seems likely. 
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The Orion Nebula (M42, NGC 1976) is th e nearest massive 



star forming region, at a distance of 450 pc ( Genzel & Stutzki 



|198£ ), thus becoming the archeytpe for studying related phe- 
nomena such as outflows, shocks and photodissociation. A 
score of planetary-mass 'bullets' were discov ered in the core 



of Ori on Molecular Cloud-One (OMC-1) by Allen fc Burton 



(1993, AB hereafter). The bullets are prominent in near-IR 



[Fell] a Fc,/2 - a -D7/2 1.64 /xm line emission, with trailing 
wakes bright in H2 1-0 S(l) 2.12 /xm line emission, which 
are sensitive to intermediate (> 50 kms~^) and slow (< 50 
kms~^) speed shocks, respectively. The bullets were iden- 
tified with fast moving optical Herbig-Haro (HH) objects 
observed by ?) in the [O I] 6300A line. As indicated by their 
line widths, they are travelling out of the cloud at speeds 
of up to 400 kms~^. Projecting back, the wakes point to an 
origin within 5 arcsec of the IRc2-complex. Their dynamical 
age is no more than 1,000 year s. A second finger system was 
found by |Stolovy et al. (199^) from HST/NICMOS obser- 



The formation mechanism for these bullet systems is 
not yet clearly understood. Two competing hypothesis are 
1) one or multiple explosive events fiing off ejecta from the 
IRc2-complex (AB) and 2) in situ Rayleigh-Tayler instabil- 
ities in time- variable winds, o r in a slow wind being over- 
taken by faster moving wind ( ^tone et al. 199*5 ). A proper 
motion study can verify the bow shock interpretation of the 
bullet-wake systems. If motions are not detected, perhaps 
because the emission arises in a standing wave in the flow 
of a wind from a central source, the ejection model (AB) 
would be invalid. Together with line profile measurements, 
proper motions may constrain the location of the central 
source by determining the three dimensional structure and 
velocities. In addition, relative proper motions of the H2 and 
[Fell] gas in a wake and bullet provide constraints on how 
ambient material is entrained by the leading bullet. 



Previously Jones & Walker (1985) measured the proper 
motion of three optical bullets (HH201, 205 and 210, fol- 



vations, within « 15 arcsec of the IRc2-complex. This inner 
flnger system emits in H2 only, with no [Fell] emission at 
their tips. These bullets are presumably interacting at lower 
speeds with the medium, leading to a lower-excitation spec- 
trum. 
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lowing the nomenclature of Reipurth 199t) by comparing 



20 plates, all taken on infrared emulsion (1-N or 103a-U) 
through an RG-8(= RG-695) filter with the Lick Observa- 
tory 120-inch telescope between 1960 and 1983. They showed 
large tangential velocities, and the authors speculated that 
the sources IRc2 and IRc9 must be responsible for the excita- 
tion, resid ing at mo st ~0.2 pc behind an obscuring molecu- 
lar cloud. Hu (1996) used narrow band [Nil 6583A images 



taken by HST/WFPC with a 4 year baseline to also de- 
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Figure 1. The [Fe II] 1.64 image of OMC-1 taken in December 1996, with proper motion vectors from the centroid method overlaid. 
For HH 117-200 (No. 12) and HH 130-119 (No. 20) the proper motion vectors from eye-estimation are used (see §3.1). The length and 
orientation of an arrow represents the transverse speed of a bullet and its direction of motion. A scale bar in the lower left corner 
corresponds to a speed of 100 kms~^ at a distance of 450 pc and its length to the distance a bullet would cover in 100 years. The 
IRc2— complex, where the driving source is presumed to be located, is at 5''35'"14!4, -5°22'30" (J2000), below the image. Every bullet 
whose proper motion has been measured is numbered (first column of Table 1). Those whose motion has S/N < 3 are marked with a 
cross only. 



termine their proper motions, and long-slit, high resolution 
echelle spectra of [O I] in measuring the three dimensional 
kinematics. He estimated shock velocities and orientation 
angles for HH201, 205, 206 and 210 adopting the bow shock 
model and analysis by Hartigan et al. (1987). He suggested 
that a hydrodynamic instability in a shell swept up by a 
poorly coUimated stellar wind from young stars is responsi- 
ble for the formation of these objects. 

Through high resol ution spectroscopic o bservations of 
the [FeH] 1.64 fj.m line, |Tedds et al. (1999a| ) obtained line 
profiles for the bullets HH 126-053 (=HH207; here HH 125- 
053, see text) and HH 120-114 (here HH 117-114). From the 
shape of the integrated [FeH] line profiles, they derived the 
bullet velocity and the orien tation angle based on the analy- 
sis of Hartigan et al. (1987 ) that showed the full- width zero 
intensity (FWZI) of the integrated profile over a single bow 
shock must be equal to the speed of the bullet itself. The re- 



sults from these studies are give in Table 1 for comparision. 
We have used images of the [Fe H] and H2 lines, taken during 
periods of exceptional seeing and with a 4-year baseline, to 
investigate proper motions of the multitude of bullets which 
are observable at near-IR wavelengths. 

Data acquisition and analysis is presented in Section 2, 
with the results from this and previous studies in Section 3. 
Discussion and plans for future work are followed in Section 
4 and 5, respectively. 



2 OBSERVATIONS & DATA ANALYSIS 

A region of ~ 2 x 2 arcmin in OMC-1 was mapped in 1992 
September 13 and 1996 November 24, givi ng a 4.2 year time 
baseline, with the IRIS near-IR camera (Allen et al. 199J) 
on the 3.9-m Anglo-Australian Telescope. Narrow-band line 
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Figure 2. The H2 1—0 S(l) 2.12 fim image of OMC— 1 taken in December 1996, with proper motion vectors overlaid as Figure |^. A newly 
found H2 bullet HH 117-212 (No. 104) is out of the frame and not shown. For its location, see Figure ^. Every H2 knot whose proper 
motions have been measured is numbered (first column of Table 1). Those with S/N < 3 are marked with a plus sign (+) only. 



images in [Fell] 1.64 /im and H2 2.12 fim emission were 
taken with 30 second integration time. A pixel scale of 0.27 
arcsec was used for both data sets with sub-arcsecond seeing 
(0.7" in 1992 and 0.9" in 1996) for both epochs. The data 
reduction process involved the usual bad pixel removal, flat- 
fielding and sky subtraction. A slight field rotation (12.2') 
between the two epochs was corrected for the second epoch 
data with a transformation solution obtained by using com- 
mon stars in the field. The error involved in image registra- 
tion is of order a tenth of a pixel. We did not attempt to 
correct for different seeing between the epochs. 

Since the bullets and trailing wakes are diffuse objects, 
their proper motion depends on whether one is measur- 
ing their bulk motion, and/or viewing changes in morphol- 
ogy as well. We adopted three methods which are sensi- 
tive to different aspects of motion; 1) image centroiding, 2) 
cross-correlation and 3) eye-estimation. Centroiding is sen- 
sitive to motion of the emission centroid of an object, and 
cross-c orrelation to morphological changes. A small area of 



of the peak intensity from the centre in this cross-correlated 
image was interpreted as proper motion. In reality, however, 
it reffects both bulk motion and changes in shape of the ob- 
ject. If the results from the first two methods do not agree, 
we attributed this to morphological changes, which was often 
confirmed by the third method - eye examination. Contours 
of the object in two epochs were compared by eye to give 
the shift in the position of intensity peaks only. 

In n aming a bullet, we fo llowed the nomenclature pro- 
posed by O'Dell & Wen (1994) to designate compact sources 
and stars in M42. It assigns to an object a catalog number in- 
dicating the position of the source. The first three digits indi- 
cate the position in right ascention (E2000), and the second 
three the position in declination. The common values for the 
inner region of M42 (5''35™ and -5°20') are not included. 
Hence, a bullet located at 5''35™11':51, -5°21'47" would 
become 115-147. The coordinates of PC043 5''35'"13°39, 
-5°21'08.0", from a HST WF/PC survey ([Prosser et al 



sky around an object is selected in two epochs and cross- 
correlated, producing a resultant image. The amount of shift 



1994 ) was used as the reference position to obtain the coor- 



dinates of bullets. 
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Figure 3. Enlargened images around HH 117-114 of H2 in 1992 {left), [Fell] in 1992 [centre] and [Fell] in 1996 {right). The compact 
[Fell] emission knot, HH 117-114 (No. 11), lies at the head of wakes in H2. Two more [Fell] bullets, HH 115-109 (No. 7) and HH 116-110 
(No. 9) are in front of it. A stationary clump 114-107 (No. 101) in H2 emission lying ahead of them. Presumably it is being rammed by 
the HH115-109/HH 117-114 (Nos. 7/11) bullet system. An [FeH] clump shows up as HH 112-109 only in 1996, and was not seen in 1992 
(middle and right panels). HH 117-057 (No. 10) is stationary, accompanyed by a bright triangle-shaped stationary H2 clump, 115-055 
(No. 102) next to it. The velocity vector and the scale bar, on the left panel, is the same as in Figure |l] H2 knots are marked by + and 
[Fe n] knots by X . 



3 RESULTS 



HH 113-156 (=HH201) & HH 117-200 [Nos. 4 & 12] 



The results of the proper motion measurements of the Orion 
bullets and wakes are listed in Table 1 (see also Figure ^ for 



[Fell] ^ nd Figure H for H2), as well as the results from pre- 



vious slLidius. The aveiage piupei muLiun uf [Ffell] buUefs 
and their trailing H2 wakes is about 8 arcsec per century for 
those whose proper motion has been measured, correspond- 
ing to a transverse speed of 170 kms~^ at the distance of 
450 pc to OMC-1. For bullets whose motion has S/N > 3, 
the average proper motion is 9.4 arcsec per century with 
a corresponding transverse speed of 200 kms~^. The direc- 
tion is similar to that inferred from their morphology. In 
general, the three methods have given similar results (aside 
from HH 109-148, HH 113-156 and HH 130-041) as there are 
few clear changes in the bullet morphology over the 4 year 
baseline. However some bullets and wakes have change in 
shape and/or intensity, as described below. We first discuss 
the proper motions seen in the [FeH] images, and then in 
the H2 images. 



3.1 Noticeable Features in [Fell] Line Emission 

Features for individual [FeH] bullets are described below in 
the ascending order of their R.A. coordinate. Their designa- 
tion from the previous studies are also given in parenthesis 
when available, and the number following is from Column 1 
in Table 1. 



This shows a double feature and/or unresolve d sub- 
structure within it 



1.5 arcsec, or ~ 600 AU in size. 



O'Dell 



& Wong (1996) identified this as four seperate objects of 



113-153, 114-155, 155-155 and 116-156. It has no counter- 
part or wake in H2, attributed to its motion close to our 
line of sight. Different methods give different proper motion 
values, mainly due to its double feature and unsymmetric 
shape. The tangential velocity Vtan calculated from centroid- 
ing, 38±21 kms~^, disagrees with those from previous stud- 
ies. This probably results from partially resolving the feature 
into two components. Vtan obtained from cross-correlation 
(177 kms~^) and from eye-estimation (100 kms~^) both 
agree reasonably. They give an orientation angle (O.A.) of 
60° for the former, and 20° for the latter method when 



combined with Hu (1996)'s estimation of the shock veloc- 
ity 310 km s^^. O.A. is the angle between the direction of 
motion of an object in space and our line-of-sight. 

HH 117-200 (No. 12) shows weak [Fell] emission, lo- 
cated just southeast of HH 113-156 (No. 4), however their 
association is not clear. The proper motion determined for 
HH 117-200 (No. 12) vary for different methods; centroid- 
ing gives a motion to the northeast with P.A.=30° while 
the other two methods to the northwest with P. A. =330°. 
Its Vtan also vatics from 90 kms~^ (eye-estimation) to 213 
kms~^ (cross-correlation). Its overall intensity distribution 
remains the same for both epochs, and a small shift in emis- 
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sion peak is attributed as the cause of this discrepancy. 

HH115-109/HH 117-114 (=HH 120-114) 

[Nos.6, 7, 9, 10 & 11] 
This bullet sys tem was designated as HH 120-114 by 
Tedds et al. (1999a ) (without exact astrometry) using AB's 
1992 image, but we seperate them into two (possibly three 
or more, see below) bullets with clear evidence of a double 
feature in the H2 wakes behind them (Figure |^) . The bul- 
lets are moving at a similar speed of ~200 kms~^ and posi- 
tion angle of 320°. HH 116-110 (No. 9), a secondary emission 
peak in the [Fell] wake of HH 115-109 (No. 7), moves slightly 
faster than HH 115-1 09 (No. 7). It is not c lear whether they 
are seperate bullets. Tedds et al. (1999a) derived the bul- 
let velocity and O.A. for this system from the shape of 
the integrated [Fell] profiles. They derived a shock veloc- 
ity Vs = 120 ± 10 kms"^ and O.A. = 60 ± 15° giving a 
corresponding Vtan of 100 ± 25 kms~^. The difference in 
vtan from that derived here indicates that using the shape 
of line profiles to infer orientation angles can be unreliable. 
In this particular case several profiles are overlapping. 

A knot at 115-103 (No. 6/103) shows weak emission in 
both [Fell] and H2 lines with a slight offset. Both species 
show similar proper motions of 100 kms~^. Note the proper 
motion in [Fell] from eye-estimation is greater than that 
from centroiding, due to a small shift in emission peak to 
west. HH 117-057 (No. 10) is found to be stationary next 
to a bright H2 clump 115-055 (No. 102). From the relative 
locations of [Fell] and H2, we suggest that this is a 'reverse' 
bullet, i.e. a dense, stationary condensation being rammed 
by outflowing material, producing a reverse-oriented bow 
shock (see §4). 

The proper motion of a knot at 118-119 was not measur- 
able because it lies in the middle of difi^use [Fe II] emission in 
the tail of the HH 115-109/HH 117-114 (Nos. 7/11) system. 
In addition the weak emission peak in 1992 had become more 
diffuse by 1996. Interestingly it lies along the axis of one of 
the following H2 wakes, suggesting this is more directly as- 
sociated with them (see §3.2). Note an [Fell] clump which 
showed up at 112-109 in 1996 (Figure |^. It was not seen in 
1992 indicating intensity variability. 

HH 116-025, HH 122-041 & HH 125-053 (=HH205, 
HH 206 & HH 207) [Nos. 5, 8, 13, 14 & 15] 

This is the most crowded field with several bullets and 
wakes overlapping (Figure Thus this field is suitable for 
determining whether any differential motion between bul- 
lets and/or between species occurs, i.e. whether or not the 
bullets move with different speeds depending on their dis- 
tance from the origin, and/or if their H2 wakes are stretch- 
ing out behind them. However, no differential motion was 
detected. This may be attributable to the short time base- 
line. If such motion was, for example, 30 kms~^, a 10-year 
baseline would be required. 

HH 116-025 (=HH205, No. 8) is one of the fastest mov- 
ing bu llets (~ 2 50kms~^) at the tip of the most crowded 
finger. Hu (1996) reported a morphological change in [Nil], 
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Figure 4. A blow-up of the most crowded field accomodating 
HH 116-025, HH 122-041 and HH 125-053 (Nos. 8, 14 & 15). The 
velocity vector and the scale bar, at the upper left corner, are the 
same as in Figure ^ 



HH 122-041 (No. 14) is irregular in shape having the same 
direction of motion with HH 116-025 and HH 125-053 (Nos. 
8 & 15). For HH 125-053 (=HH207, HH 126-053; No. 15), 
~ 100 k ms~^ proper motion is measured while the estima- 
tion by Tedds et al. (1999c) gives vtan = 140 ± 25kms~^, 
with Vs = 150 ± lOkms"' and O.A. = 70 ± 15°. Emission 
at the bullet head is observed to get stronger without any 
morphological change. 



but no change is apparent in our [Fe II] image. His estimation 
of Vs ~ 310kms~^ with O.A. — 90° (moving on the plane- 
of-sky), makes the expected tangential velocity the same as 
the shock velocity. A smaller bullet, HH 114-023 (No. 5), lo- 
cated further northwest of it, shows similar proper motion. 



HH 130-119 [No. 20] 

HH 130-119 is located along the line of a series of promi- 
nent bullets; HH 116-025, 122-041 and 125-053 (Nos. 8, 14 
and 15, respectively). It is quite diffuse in shape with no 
noticeble emission peak in the 1992 image. In 1996, how- 
ever, an emission peak was prominent in the lower-left cor- 
ner of the bullet, considerably affecting the proper motion 
measurement using image centroid and cross-correlation. By 
blinking between two epoch data, we can see its motion to- 
ward northwest (P.A.=5°), and this value is used in Figure 1. 
This demonstrates the limitation of centroiding when a bul- 
let changes its intensity distribution over time. Apparently 
it has the largest proper motion, but this is a reflection of 
the change in distribution within its diffuse, irregular mor- 
phology, rather than of actual bulk motion. 

HH 154-041 (=HH 210) [No. 27] 

Jones & Walker (lyso) reported a 'nucleus' in its head. 
We also see this in [Fell] and in J band (1.25 /im), as much 
fainter emission with similar morphology, but not in H2 
emission. This could mean that its motion is close to our 
line of sight, or that it has escaped from the molecular cloud, 
hence no H2. While zero proper motion was expected from 
Hu (1996| )'s estimation of Vs = 400 kms'^ and O.A. = 0° 
(moving directly toward us), our measurement gives vtan ~ 
260 kms"'^. In fact, the morphology suggests there are at 
least two bow features overlapping - while the overall shape 
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Figure 5. HH 154-041 (=HH210, No. 27) is shown in [Fell] line 
emission with the proper motion vector overlaid. Two or more 
bow features overlap here. While the overall shape is preserved, 
the lower emission peak has slightly shifted to the right, in com- 
parison to the motion of the bullet head. The velocity vector and 
the scale bar, on right panel, are the same as in Figure 1. 



[Fell] Bullets 




-05 22 12 

5''35'"13'0 



12r5 12?0 1ir5 
RA (J2000) 



1 I'O 5"35'"13'0 



12S 12M lirs 

RA (J2000) 



Figure 6. (left) Newly discovered H2 bullets arc shown with 
proper motion vectors overlaid as Figure 1. They are found be- 
tween HH 113-156 (=HH201, No. 4; right) and the presumed ori- 
gin in the IRc2-complex. The corresponding [Fe II] line emission 
is shown right. The proper motion vector for HH 113-156 (No. 4) 
is offset for clarity, and the proper motion vector from the eye- 
estimation is used for HH 117-200 (No. 12; see text). The velocity 
vector and the scale bar, on right panel, are the same as in Figure 
1. 



is preserved the second emission peak shows a small shift 
(0.6") to southwest, different from the northeast motion of 
the tip of the bullet. 



3.2 Discovery of New H2 Bullets 

We report the discovery of new H2 bullets; HH 117-212, 
HH 119-203 and HH 120-145 (Nos. 104, 105 and 107, re- 
spectively; Figure ^). They emit in H2 only, located be- 
tween HH 113-156 (=HH201, No. 4) and the presumed origin 



in the IRc2-complex. Their morphology resembles that of 
[Fell] bullets but with smaller proper motions of order '^lOO 
kms~^. We presume their slower speed results in a lower- 
excitation spectrum without [FeH]. These are thus similar 
to the inner system of H2 bullets discovered within 15 arcsec 



of IRc2-complex by Btolovy et al. (1998 



3.3 Noticeable Features in H2 Line Emission 

The proper motion measurements for both H2 emitting 
wakes and H2 bullets were difficult to determine because 
of their diffuse morphology. H2 wakes which do have leading 
[Fe H] bullets appear to have a similar proper motion vector 
as their leading bullet. However, in general, the proper mo- 
tions of H2 bullets are slower than those of [FeH] bullets. 

HH 114-107, 115-055 & 115-103 [Nos. 101, 102 & 103] 

HH 114-107 (No. 101) at the tip of HH 115-109/117-114 
(Nos. 7/11) has no associated [Fell] emission (Figure ^) and 
is found to be stationary. This may be a stationary cloudlet 
being rammed by the shock wave running ahead of the HH 
obje ct, as is the case with t he 'tadpoles' in Helix Nebula 
(eg. Burkert & O'Dell 1998). In such a case, the nearside 



of the clump to the IRc2-complex would show the highest 
excitation with its trail pointing away from the complex. 
Unfortunately the size of the knot is not big enough for 
such structure to be resolved. 

The same explanation can be applied to 115-055 
(No. 102) which is a bright triangle-shaped clump with a 
weak [Fell] emission knot HH 117-057 (No. 10) interior to it. 
When their emission peaks in [Fe II] and H2 are connected, 
the line has P. A. of ~ 310°, pointing about 30 arcsec east of 
the IRc2-complex. 

For 115-103 (No. 103), the emission peak has shifted to 
west by half an arcsecond giving much larger proper mo- 
tion from eye-examination than from centroiding, as has its 
counterpart in [Fell] emission HH 115-103 (No. 6). 

120-120 [No. 106] 

Two bow features are discernable in H2 wakes associ- 
ated with an emission peak at 120-120 (No. 106) (Figure 
0). They are trailing the [Fell] bullets HH 115-109 (No. 7) 
with P.A.=330° and HH 117-114 (No. 11) with RA.=320°. 
A line, drawn along the presumed axis of the smaller bow 
feature, passes through a secondary [Fell] emission peak 
at 118-119 which coincides with the tip of the bow. Ex- 
tended further northwest, the line passes through HH117- 
114 (No. 11). These features may arise from the fragmenta- 
tion of the bullet into two pieces, rather than these being 
two seperate bullets projected on top of one another. 

Figure ^ also shows a morphological change at the head 
of these wakes. The cone-shaped tip of the H2 wake in 1992 
has disappeared leaving it flat-headed by 1996. Inhomogene- 
ity in the ambient molecular cloud on scales of 100 AU, 
causing changes in the rate of entrainment of material be- 
hind a bullet, could produce such a change. 



4 DISCUSSION 

The most basic result from the data presented here is that 
the bullet/bow-shock interpretation placed on the observa- 



An infrared proper motion study of the Orion bullets 7 




5"35"'12!'0 11=8 11=6 11=4 
RA (J2000) 



5''35'^12=0 11=8 11=6 11=4 
RA (J2000) 



Figure 7. The open head of wakes accomodating H2 clump 
120-120 (No. 106), lead by the [Fell] bullet system HH115- 
109/HH 117-114 (Nos. 7/11). A straight line is drawn along the 
presumed axis of the smaller bow feature with P. A. of 320° . A star 
sign marks the position of a [Fell] knot 118-119, and X the [Fell] 
bullets HH 115-109 (No. 7), HH 116-110 (No. 9) and HH 117-114 
(No. 11). The velocity vector and the scale bar in the right panel 
are the same as in Figure 1. See Figure ^ for the proper motion 
vectors for the [Fe II] bullets shown here. 



tions of the [Fell] caps and H2 wakes in Orion by AB is 
essentially correct. The proper motions measured confirm 
that the bullets are travelling away from the core of the 
cloud at supersonic speeds, and are thus projectiles that are 
being shocked. They cannot be stationary cloudlets, nor can 
they be standing waves in a flow. 

Their speeds, typically of order 200 kms~^, and di- 
rections are consistent with both the widths of their line 
profiles, and with the morphology of individual bullets and 
fingers. We are observing bullets and associated wakes inter- 
acting with an ambient medium, spread over a wide opening 
angle, as they move outwards from the centre of the molec- 
ular cloud. 

The inferred speeds are much in excess of the dissocia- 
tion speeds for H2 and even [Fell], and thus cannot be the 
actual shock speed of the line emitting gas. There is, how - 
ever, no suggestion from other data (eg. Tedds et al. 1999h) 
that the bullets are traversing a moving medium (thus re- 
ducing the shock velocity to the difference between the bul- 
let and medium velocities). Oblique shocks, except at the 
very head, must occur with a significant tangential compo- 
nent of velocity remaining in the shocked material. This is 
as expected for a bow shock around a bullet, with only the 
normal component of the bullet's motion at each point along 
the bow contributing to the shock speed. 

We see little evidence for significant morphological 
changes in the bullets over the four-year baseline, nor for 
fading or brightening of individual knots over the same 
time, aside from [Fell] clumps HH 112-109 and HH 130- 
119 (No. 20), and the head of the H2 wake associated with 
120-120 (No. 106) (see §3). This is thus different from the 
observation o f some of the bull ets in HH46/47 (?) and 
HH 111/121 (Coppin et al. 1998), where brightness varia- 



tions on timescales of 1 year (the cooling time for H2) are 
seen. Such variations likely arise from inhomogeneities in 
the ambient medium over spatial scales of at most some 
tens of AU, leading to changes in the amount of material 
being swept up over this time. For the bullets in Orion the 
medium must thus be reasonably homogeneous over ~ 100 



AU scales, otherwise changes in shock strength, and thus 
line intensities, would be observed. However, the change in 
H2 wake intensity near 120-120 (No. 106) and appearance of 
[FeH] clump HH 112-109 (see Figure ^) suggests this region 
is not homogeneous over this scale. Here, a bullet may have 
fractured into several pieces, each forming its own set of H2 
wakes, overlapping behind. The similar transverse speeds 
of [FeH] bullets HH 115-109, HH 116-110 and HH 117-114 
(Nos. 7, 9 & 11) support this suggestion. 

No differential motion, however, between [Fell] bullets 
and their H2 wakes is observed. While limits on such motion 
remain weak over a four year timescale, this suggests that 
the bullets are not separating from their wakes. Together 
with the lack of any morphological changes, the wakes also 
cannot be stretching. In other words, each bullet and wake 
appears to be moving together. This does not, however, im- 
ply that the gas within the bullet /wake is stationary with 
respect to the bullet velocity. Material swept-up by the bul- 
let will be swept down the wake at the tangential compo- 
nent of the shock velocity, and will thus traverse the length 
of the bow. It will continue to be observed until it reaches 
the Mach point, where the normal component of the bullet 
velocity drops below the minimum shock speed to excite the 
gas (~ 5 kms~^ for molecular gas). The time taken for this 
will be of order 100 years. Since the wakes are not stretching 
either, the bullets must also have been in motion for at least 
this long, in order to reach a steady-state configuration. 

The bullets observed in H2 only, on the inner edge 
of the finger system HH 113-156 (No. 4), are moving more 
slowly 100 kms~^) than those emitting in [Fell] (~ 200 
kms~^), though their speeds still exceed that for dissocia- 
tion of H2. However, these slower speeds are consistent with 
lower excitation for the inner bullets, as would be inferred 
by the absence of [Fell] in them. This gives support to the 
inner H2 bullet system discovered by Stolovy et al. (1998) 
within 15 arcsec of IRc2-complex as being slower moving 
analogues of the [Fell] bullets. 

Nevertheless not all line-emitting knots are moving out- 
wards at these speeds. For instance, HH 117-057 (No. 10) and 
115-055 (No. 102) appears to be stationary with respect to 
the ambient cloud. HH 114-107 (No. 101) is also stationary 
with HH115-109/HH 117-114 (Nos. 7/11) moving towards 
it. We hypothesise that it is a dense clump in the ambi- 
ent medium just now being rammed by HH 115-109/HH 117- 
114 (Nos. 7/11). Time variability is evident in [Fell] clump, 
HH 112-109 in Figure ^, which showed up in 1996 but was 
not seen in 1992. In all these cases, we expect the clumps to 
brighten over the next few years as they are further overrun. 

The proper motion vectors of the bullets clearly extend 
back towards an origin somewhere in the vicinity of IRc2- 
complex. However their extrapolation remains no more ac- 
curate than that undertaken by Allen & Burton (1993), who 
extended the wakes back towards an origin within 5 arcsec of 
the IRc2-complex. Our present data still cannot shed any 
more light on where this might be, which awaits a longer 
time baseline of measurements. 

The projected distance of bullets from IRc2-complex 
is plotted against its proper motion speed in Figure |^. Ex- 
cluding HH 130-119 (No. 20) which was problematic in its 
measurement, the plot suggests that more distant bullets, 
in general, are moving faster. This is consistent with an ex- 
plosive, or at least an impulsive event, as the origin for the 
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0.1 0.2 0.3 0.4 

Distance [Parsec] 

Figure 8. The bullet speeds whose S/N > 3 are plotted against 
their projected distance from the IRc2-complex, again numbered 
with those in the first column of Table 1 . Thick solid lines are for 
HH 116-025 (No. 8), 122-041 (No. 14) and 125-053 (No. 15) from 
right to left (see §4). 



bullets. Certainly it must have occurred over a timescale 
much less than that of the expansion, i.e. <^ 1,000 years. 
This is most clearly demonstrated for HH 125-053, 122-041 
and 116-025 (No. 8, 14 and 15, respectively) along the main 
finger, which show a linear increase in speed with distance. 
However the motions are not consistent with free expansion 
from the IRc2-complex; extrapolation to zero speed would 
suggest an origin 0.05 - 0.1 pc from the IRc2-complex. This 
might correspond to where the bullets have formed, through 
instabilities in a time variable wind. 



5 FUTURE WORK 

Our data represent just two epochs of observation, and is 
only the start of an investigation of the motion of the Orion 
bullets. With further epochs of data we will be able to ad- 
dress further questions regarding their origin and excitation, 
including: 

(a) whether the bullets are decelerating (changing in 
velocity)? 

(b) are there non-radial motions, perhaps imparted by 
collisions with dense objects? 

(c) are there more differential motions between the bul- 
lets and their wakes? 

(d) could we be seeing a pattern speed of a bullet, rather 
than its bulk motion? This may be evident if wiggles are 
observed on a proper motion vector. 

(e) are there more morphological or intensity changes in 
time, such as the stretching of a wake or the brightening or 
fading of a knot? These might arise from inhomogeneities in 
the ambient medium, thus indicating what scale these occur 
on. 

Further proper motion vectors will help tie down the 
origin of the driving source, by reducing the errors on their 
extrapolation back to an origin. Combining the data with 
line profile measurements would give the 3D velocity of a 



bullet, and allow a 3D picture of their distribution in space 
to be built up. The line profiles could then b e compared to 
models such as those of Hartigan et al. (1987) which predict 
profile shapes for a range of orientation angles. In addition, 
determining an orientation directly will help constrain shock 
models for the line production, such as the form of the cool- 
ing law and shape of the bow, which significantly effect the 
resultant profile. The [FeH] data can also be compared to 
optical data obtained with HST, which should yield more 
accurate velocities on shorter time base lines due to the su- 
perior spatial resolution that can be obtained. 
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Table 1 . Proper Motions of the Orion Bullets and Wakes in [Fe II] and H2 emission. 
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The proper motion of bullets in [Fell] is followed by the proper motion of bullets and wakes in H2 (from No. 101 to 111). In naming H2 
knots, HH indicating Herbig-Haro object is only given when it is an H2 bullet, and not when it is the wake behind an [Fell] bullet. 

^ Prom the shift in the intensity centroids of a bullet in two epochs. A proper motion of 10 arcsec per century corresponds to a transverse 
speed of 215 kms~^ at the distance of 450 pc to OMC— 1. 

2 The tangential velocity; the velocity projected onto the plane-of-sky. 

^ The position angle; the direction of motion in the plane-of-sky, starting from N and increasing to E. 

* From the cross— correlation of the distribution of emission from the bullets at two epochs. It generally agrees with the results from the 
centroid (see §2). The diffuse shape of H2 knots sometimes makes cross— correlation unreliable, when their value is not given. 

^ From eye— estimation of the contour maps of bullets. Eye-estimation is sensitive to changes in the intensity peak only. 



^ The s hock velocity and orientation angle. The shock velocities are estimated from the line profile of [OI] (Hu 1996) and [Fell] (Tedds 



et al. IE 99a). The orientation angle is the angle between our line-of-sight and the direction of motion of an object. O.A.=90° when the 
object is moving on the plane-of-sky. 

* For these bullets, the proper motion vectors from the eye-estimation (P.A.=333° for No. 12 and 5° for No. 20, respectively) are used in 
Figure ^ (see text). 

* See §3.2 for the discrepancy in the proper motion values from different methods. 

(a) IHu (19961), d=440 pc used; (b) Ijones & Walker (19851), d=460 pc used; (c) lledds et al. (1999a|) 



